MAGNETIC PROPERTIES OF CUPRATE 
PEROVSKITES 



A. Sherman 

Institute of Physics, University of Tartu, Riia 142, 51014 Tartu, Estonia 
alexei@fi.tartu.ee 



Abstract The magnetic susceptibility of underdoped yttrium and lanthanum cup- 
rates is interpreted based on the self-consistent solution of the t-J model 
of a Cu-O plane. The calculations reproduce correctly the frequency de- 
pendencies of the susceptibility in YBa2Cu307_ y and L^-^Sr^CuCU 
attributing their dissimilarity to the difference in the damping of spin 
excitations. In YBa2Cu307_ !/ these excitations are well defined at the 
antiferromagnetic wave vector Q = (rr, tt) even in the normal state 
which manifests itself in a pronounced maximum - the resonance peak 
- in the susceptibility. In L^-xSr^CuCU the spin excitations are over- 
damped which leads to a broad low-frequency feature in the suscepti- 
bility. The low-frequency incommensurability in the magnetic response 
is attributed to a dip in the magnon damping at Q. The calculated 
concentration and temperature dependencies of the incommensurabil- 
ity parameter conform with experimental observations. Generally the 
incommensurate magnetic response is not accompanied with an inho- 
mogeneity of the carrier density. 
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1. Introduction 

Among the results obtained with inelastic neutron scattering is the 
detailed information on the magnetic susceptibility in YE^C^C^-,, 
and lanthanum cuprates which reveals a considerable difference in their 
magnetic properties. A sharp magnetic collective mode - the so-called 
resonance peak - was observed at frequencies u r = 20 — 40 meV in 
YBa2Cu307_ y and some other cuprates, in underdoped case both be- 
low and above T c [1]. In the momentum space the mode is strongly 
peaked at the antiferromagnetic wave vector Q = (tt, tt). Contrastingly, 
no resonance peak was observed in lanthanum cuprates. Instead for 
low temperatures a broad feature was detected at ujf ~ 7 meV and for 
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this and lower frequencies the susceptibility was found to be peaked at 
incommensurate momenta (tt ± 2ir5, it) and (tt,tt ± 2tt5) [2]. Recently 
the low-frequency incommensurate magnetic response was also observed 
in YBa2Cu307_ y [3]. Presently the fundamental difference of the sus- 
ceptibility frequency dependencies in these two groups of cuprates is 
actively debated, as well as the nature of the resonance peak and the 
low-frequency magnetic incommensurability. 

The aim of this work is to demonstrate that the above-mentioned 
unusual properties of cuprates can be interpreted in the framework of the 
t-J model of a Cu-0 plane which is a common structure element of these 
crystals. The model was shown to describe correctly the low-energy part 
of the spectrum of the realistic extended Hubbard model [4]. To take 
proper account of strong electron correlations inherent in moderately 
doped cuprate perovskites the description in terms of Hubbard operators 
and Mori's projection operator technique [5] are used. The self-energy 
equations for hole and spin Green's functions obtained in this approach 
are self-consistently solved for the ranges of hole concentrations < x < 
0.16 and temperatures 2 K< T <1200 K. Lattices with 20x20 sites and 
larger are used. 

The calculations reproduce correctly the frequency and momentum 
dependencies of the resonance peak in YBa2Cu307_ J/ , and its variation 
with doping and temperature in normal and superconducting states. 
The peak is connected with the excitation branch of localized Cu spins 
and the peak frequency is close to the frequency of these excitations at 
Q. The absence of the resonance peak in lanthanum cuprates is related 
to an increased damping of spin excitations which is possibly connected 
with a large hole damping in these crystals. For low frequencies and 
temperatures the susceptibility is peaked at incommensurate wave vec- 
tors (-7T ± 2tt5, 7t) and (w, ir ± 2n5). The incommensurability is connected 
with a dip in the magnon damping at Q. In agreement with experiment 
the incommensurability parameter 5 is approximately proportional to x 
for 0.02 < x < 0.12 and saturates for larger concentrations. The incom- 
mensurability disappears with increasing temperature. Generally the 
incommensurate magnetic response is not accompanied by an inhomo- 
geneity of the carrier density. 

2. Main formulas 

The Hamiltonian of the two-dimensional t-J model reads 



H — inm ffl L fl mff + „ Jnm ( s n s m + s n s m ) > (-Q 



Magnetic Properties of Cuprate Perovskites 



3 



where a n<J = |ncr)(nO| is the hole annihilation operator, n and m label 
sites of the square lattice, a = ±1 is the spin projection, \na) and 
|n0) are site states corresponding to the absence and presence of a hole 
on the site. These states may be considered as linear combinations 
of the products of the 3d x 2_ y 2 copper and 2p a oxygen orbitals of the 
extended Hubbard model [6]. We take into account nearest neighbor 
interactions only, t nm = -t^ a ^ n ,m+a and J nm = JJ2 a 5 n ,m+a where 
the four vectors a connect nearest neighbor sites. The spin-^ operators 
can be written as s„ = \ J2a c|ncr)(ncr| and s£ = |n<j)(n, — <r\. 

Properties of the model are determined from the hole and spin re- 
tarded Green's functions 

G(kt) = -i9(t)({*Ut),<b})> = -i*(*)K(t),*ik]>, (2) 

where Oko- and s k are the Fourier transforms of the respective site op- 
erators, operator time dependencies and averaging are defined with the 
Hamiltonian TC = H — /i ^ n ai„a n(J with the chemical potential fi. As 
mentioned, to obtain self-energy equations for these functions we used 
Mori's projection operator technique [5, 7]. In this approach the Fourier 
transform of Green's function {(Aq\Aq)) is represented by the continued 
fraction 

«4>l4» = • ( 3 ) 

LO - E y- 

uo — E\ 



The elements of the fraction E{ and V{ are determined from the recursive 
procedure 

[A n , H\ = E n A n + A n+ \ + Vn-iAn-i, 

E n = \[A n ,H]-Al\\A n -Ai\-\ (4) 
V n - 1 = \A n -Ai\\A n ^ 1 -Al_ 1 \- 1 , V-! = 0, n = 0,1,2,... 

The operators A, constructed in the course of this procedure form an or- 
thogonal set, | At -At | cx Sij. In Eqs. (3) and (4) the definition of the inner 

product | Ai ■ A'j | depends on the type of the considered Green's func- 
tion. For example, for functions (2) these are ({Ai,A^}) and (LAj, Aj]), 
respectively. The method described by Eqs. (3) and (4) can be straight- 
forwardly generalized to the case of many-component operators which is 
necessary, for example, to consider Green's functions for Nambu spinors 
in the superconducting state [7]. 



The residual term of fraction (3) is the Fourier transform of the quan- 
tity 

T=\A nt -Ai\\A n _ 1 -A\ l _ l \-\ (5) 

where the time evolution of the operator A n is determined by the equa- 
tion 

d ™ _ i 

l J t Ant = II^ 1 - Pk)[A n t,H], A nyt=0 = A n (6) 
/• o 

with the projection operators P n defined as P n Q = \Q-A' n \ \A n -Al l \~ 1 A n . 
The residual term T is a many-particle Green's function which can be 
estimated by the decoupling. The decoupling is also used in calculating 
the elements of the continued fractions (3). Following Ref. [8] this ap- 
proximation is improved by introducing the vertex correction a which is 
determined from the constraint of zero site magnetization 

(si) = \(l-x)-(s n 1 st 1 ) = 0. (7) 

Since the statistical averaging includes samples with different ordering 
orientations, this condition is fulfilled also in magnetically ordered states. 

Equations for Green's functions obtained in this way for the normal 
state read [7] 

GN = wT~^' D kw = 2 ^7T~\ 2~' 

16vrf 2 r°° 



1 + 7k' 
1 - 7k' 



oo 

2 



7k-k' + 7k + sgn(w')(7k-k' - 7k) 
[riB(-w') + n F (u) - oj')} 



Imn(ko;) 



(8) 

x A(\l-\s!,u-J)B(\s!J), 

97Tt 2 J 2 (l - x) V"^, , 2 

I -.Vfl^UfEl 2^(lk+V - 7k') 



2N(l-j k )(J\C 1 \+tF 1 ) 

/oo 
dw'A(kV)A(k + k',ci; + w / ) 
-oo 



k' 

, n F (oj + uj') - n F {u)') 



UJ 



where <f> = \{1 + x), 7 k = ±[cos(k x ) + cos(k y )], C\ = (s„ ls n+a) an ^ 
F\ = (on a n+a) ar e the spin and hole correlations on neighboring sites 
which, as well as the hole concentration x, can be calculated from Green's 
functions, 

col = 16J 2 a|C 1 |(l- 7 k)(A+l+7k), e k = -(4^+6df H+ZFrf- 1 J) 7k 

(9) 
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are the energies of spin excitations and holes with the parameter A 
describing the spin gap at Q, A(kw) = — 7r _1 Im G(kw) and B(\ao) = 
— 7T _1 Im D(kw) are the hole and spin spectral functions, N is the number 
of sites, np{ui) = [exp(cj/T) + l] -1 , ng(w) = [exp(cj/T) — l]" 1 with 
the temperature T. Real parts of the self-energies are calculated from 
their imaginary parts and the Kramers-Kronig relation. The self-energy 
equations for the superconducting state can be found in Ref. [7]. 

Equations (7)-(9) form a closed set which can be solved by iteration 
for fixed chemical potential \x and temperature T. We carried out such 
calculations for the parameters of the model J = 0.1 eV, t = 0.5 eV 
which correspond to hole-doped cuprates [9]. To stabilize the iteration 
procedure an artificial damping irf, r] = (0.015 — 0.05)t, was added to 
the frequency in the hole Green's function. 

The results of such calculations can be directly compared with the 
imaginary part of the magnetic susceptibility x" derived from the neu- 
tron scattering experiments, since this quantity is connected with B(kw), 

»/, x 2 „/, n 16/4(1 -7k)(«/|Ci| +tF 1 )ujlmU . . 

where (1b is the Bohr magneton. Equation (10) is the most general 
form for the magnetic susceptibility which follows from the Mori formal- 
ism. With the spin excitation damping II(ku;) described by the fermion 
bubble (8) this equation acquires some similarity with the susceptibility 
derived in the random phase approximation [10] which is frequently used 
for the discussion of the magnetic properties of cuprates. For this similar- 
ity the term in the denominator of Eq. (10) has to be negligibly small 
in comparison with u;Ren(ku;). For the hole spectrum self-consistently 
calculated here this can happen only in the overdoped region x > 0.3 
when the spin correlation C\ which determines the magnitude of ui^ [see 
Eq. (9)] tends to zero [7, 11]. 

3. The spectrum of spin excitations 

As follows from Eq. (10), the frequencies of spin excitations satisfy 
the equation 

u? - wRen(kw) -ul = 0. (11) 

As seen from Eq. (8), Im II(ka;) is finite for k — > 0, whereas Wk vanishes in 
this limit. Therefore the spin Green's function has a purely imaginary, 
diffusive pole near the T point, in compliance with the result of the 
hydrodynamic theory [12]. 

The dispersion of spin excitations near Q is of special interest, because 
this region gives the main contribution into the neutron scattering. This 
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dispersion is shown in Fig. 1. As seen from this figure, in contrast to the 




Figure 1. The dispersion of spin excitations near Q along the edge of the Brillouin 
zone. The dispersion was calculated in a 20x20 lattice for x = 0.06 and T = 17 K 
(filled squares, the solid line is the fit with oik = [<^q + c 2 (k — Q) 2 ] 1 ^ 2 ). Open squares 
are the experimental dispersion [1] of the maximum in the frequency dependence of 
the odd x"(qw), q = k - Q, in YBa 2 Cu 3 6 . 5 (a; « 0.075 [13]) at T = 5 K. 



usual spin-wave theory the frequency of spin excitations at the antiferro- 
magnetic wave vector Q is finite. This frequency ojq = 4J(2a|Ci (A) 1 / 2 
is directly connected with the spin correlation length £. Indeed, using 
Eq. (8) and taking into account that the region near Q gives the main 
contribution to the summation over k, we find for large distances and 
low temperatures 

(«> = N- 1 £e ikl f^du coth (^) B(kw) oc e^\i/\\\) 1/2 ^ m , 

(12) 

where £ = ^aA 1//2 and a is the distance between in-plane Cu sites. As 
follows from the calculations, for low concentrations and temperatures 
A « 0.2a; and consequently £ « ax -1 / 2 . An analogous relation between 
the spin correlation length and the concentration has been derived from 
experimental data in La2_ x Sr :r Cu04 [14]. 

In accord with the Hohenberg-Mermin- Wagner theorem [15] the con- 
sidered two-dimensional system is in the paramagnetic state for T > 0. 
This result can also be obtained using Eq. (8). Also it can be shown 
[7] that in the infinite crystal at T = the long-range antiferromagnetic 
order is destroyed for the hole concentrations x > x c w 0.02. 
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4. The resonance peak 

The imaginary part of the spin susceptibility at the antiferromag- 
netic wave vector, obtained in our calculations, is shown in Fig. 2. The 
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Figure 2. The imaginary part of the spin susceptibility at Q in the ri-wave super- 
conducting (a,b) and normal (c,d) states. Curves show the results of our calculations 
for r] = 0.015i, T = 17 K (a,b) and T = 116 K (c,d), x = 0.06 (a,c) and x = 0.12 
(b,d). Squares are the odd susceptibility measured [1] in YBa2Cu30e.5 (T c = 45 K, 
x « 0.075, a and c) and in YBa 2 Cu 3 6 .83 (T c = 85 K, x « 0.14 [13], b and d) at 
T — 5 K (a,b) and at T = 100 K (c,d). To depict the experimental data in the same 
scale they were multiplied by the factor 1.5. 

experimental data on \" m YBa2Cu307_ 2/ [1] are also depicted here. 
YBa2Cu307- 3/ is a bilayer crystal and the symmetry allows one to di- 
vide the susceptibility into odd and even parts. The odd part, which for 
the antiferromagnetic intrabilayer coupling can be compared with our 
single-layer calculations, is shown in the figure. 

The maximum in Fig. 2 is the resonance peak. As seen from this 
figure, the model reproduces satisfactorily the evolution of the peak po- 
sition and shape with doping and temperature. In the t-J model the 
peak is connected with the coherent excitations of localized Cu spins 
near the antiferromagnetic wave vector Q and the peak frequency for 
k = Q is close to the frequency wq in Fig. 1. It is also seen in this figure 
that the experimental dispersion of the peak is close to the dispersion of 
the mentioned spin excitations. Thus, we came to conclusion that the 
resonance peak in YBa2Cu307_, / and apparently in other cuprate per- 
ovskites where it was observed is a manifestation of the magnon branch 
modified in the short-range antiferromagnetic order. 
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The above results are related to the underdoped case where the reso- 
nance peak is observed both in the normal and superconducting states. 
In the overdoped case in YBa2Cu3C , 7_j / the peak is observed in the su- 
perconducting state only [1]. Calculations for the t-J model show that 
in the normal state on approaching the overdoped region the maximum 
in x" rapidly loses its intensity, broadens and shifts to higher frequencies 
[16]. This result correlates with the broad feature observed experimen- 
tally in these conditions. The mechanism of the peak reappearance in 
the superconducting state in the overdoped region was considered in 
Ref. [17]. The opening of the d-wave superconducting gap with the 
magnitude 2A S > loq decreases considerably the spin excitation damp- 
ing near Q which restores the peak in \" ■ 

As seen from Fig. 2, the resonance peak has a low-frequency shoulder 
which is more pronounced for low concentrations and temperatures. The 
shoulder stems from the frequency dependence of the magnon damping 
ImII(Qa;) which in its turn is a consequence of the hole Fermi surface 
nesting existing in the t-J model for moderate doping [7]. An analogous 
nesting is expected in the two-layer YBa2Cu307__ J/ between the parts of 
the Fermi surface belonging to the bonding and antibonding bends [10]. 

It is worth noting that for all four curves in Fig. 2 the value of 
|II(Q,u;q)|/2 is smaller than cjq. The satisfactory agreement between 
the experimental and calculated results allows us to conclude that near Q 
the spin excitations are not overdamped in underdoped YBa2Cu3C , 7_j / . 
In calculating the curves in Fig. 2 the artificial broadening r] in the 
hole Green's function was set equal to 0.015t. This broadening simu- 
lates contributions to the hole damping from mechanisms differing from 
the spin excitation scattering. The magnitude of these interactions can 
essentially vary in different crystals. The example of the magnetic sus- 
ceptibility calculated with the larger hole damping rj = 0.05i is shown in 
Fig. 3. As seen from this figure, the spin excitation damping is very sen- 
sitive to the hole damping. The increase of this latter damping leads to a 
marked growth of the spin excitation damping which results in the over- 
damping of these excitations. Instead of the pronounced peak at the fre- 
quency ujq a broad low-frequency feature is observed in the susceptibil- 
ity (see Fig. 3b). This shape resembles that observed in L^-xSr^CuO^ 
[2]. The comparison was carried out in Fig. 3a. It should be noted 
that the use of a comparatively small 20x20 lattice did not allow us 
to describe the incommensurability of the magnetic response - x" is 
peaked at Q and our calculated data belong to this momentum (the use 
of a larger lattice is too time-consuming for the self-consistent calcula- 
tions). In Lai.86Sro.i4Cu04 the low-frequency susceptibility is peaked 
at incommensurate momenta k = (it ± 27r<5, it), (it, it ± 27r<5) and the ex- 
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Figure 3. The frequency dependence of the imaginary part of the susceptibility, (a) 
The solid line corresponds to T = 29 K, r\ — 0.05i and x « 0.1. The dashed line is for 
T = 116 K and the same other parameters. Symbols are experimental results [2] in 
Lai.86Sro.i4Cu04 at T = 35 K (filled circles) and T = 80 K (open circles). Vertical 
bars show experimental errors. The calculated susceptibility is given for k = Q, the 
experimental data are for a wave vector of the incommensurate peak, (b) The solid 
line is the same as in part (a). The dashed line is calculated for n — 0.015t, other 
parameters are the same as for the solid line. 



perimental data [2] shown in Fig. 3a belong to one of these momenta. 
These data were increased by approximately 2.5 times to show them in 
the scale with the calculated results. The need for scaling is apparently 
connected with the splitting of the commensurate peak into the four 
incommensurate maxima. Apart from this difference in the momentum 
dependencies, the calculated frequency and temperature dependencies 
are in good agreement with the experimental results. Thus, it can be 
concluded that the dissimilarity of the frequency dependencies of the sus- 
ceptibility in YBa2Cu307_ ?/ and La2- x Sr x Cu04 is connected with differ- 
ent values of the damping of spin excitations, which are well defined at 
the antiferromagnetic wave vector in the former crystal and overdamped 
in the latter. This property of spin excitations in L^-zSr^CuC^ is not 
changed in the superconducting state due to the small superconducting 
gap 2A S ss 9meV < ujq in this crystal [17]. 

5. Incommensurability in the magnetic response 

For low frequencies susceptibility (10) is essentially simplified, 

„ u;Imn(ku;) 

x (M °c 4 — • ( 13 

As seen from Fig. 1 and Eq. (9), ui^ is the increasing function of the dif- 
ference k — Q and therefore the denominator of the fraction (13) favors 
the appearance of a commensurate peak centered at the antiferromag- 
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Figure 4- The Brillouin zone of the square lattice. Solid curves are two of four 
ellipses forming the Fermi surface at small x. Dashed lines are the Fermi surface 
contours shifted by ±(k — Q). Regions of k' and k + k' contributing to the damping 
(14) are shaded. 



netic momentum in the susceptibility. However, if the spin excitation 
damping has a pronounced dip at k = Q the peak splits into several 
peaks shifted from Q. To make sure that the damping may really have 
such a behavior let us consider the case of low hole concentrations and 
temperatures. In this case the hole Fermi surface consists of four ellipses 
centered at (±7r/2, ±vr/2) [7]. Two of them are shown in Fig. 4. The spin 
excitation damping described by the fermion bubble (8) is simplified to 

ImII(ku;) oc ^(7k' - 7k'+k) 2 »H £ k'+k - A*) ~ «iK £ k' ~ m) 

k' 

x S(lo + e k / - e k '+k) (14) 

in the considered case. Here £k is the hole dispersion. In the process 
described by Eq. (14) a fermion picks up an energy and momentum of 
a defunct spin excitation and is transferred from a region below the 
Fermi level to a region above it. However, for k = Q this process is 
impossible, because for initial states k' interior to an ellipse in Fig. 4 all 
final states with momenta k + k' will be inside another ellipse, i.e. below 
the Fermi level. Thus, in this simplified picture the damping vanishes 
for k = Q and grows with increasing the difference k — Q approximately 
proportional to the shaded area in Fig. 4 [18] . 

With increasing the hole concentration the Fermi surface of the t- 
J model is transformed to a rhombus centered at Q [16]. This result 
is in agreement with the Fermi surface observed in La2- x Sr x Cu04 [19] 
[however, to reproduce the experimental Fermi surface terms describing 
the hole transfer to more distant coordination shells have to be taken 
into account in the kinetic term of Hamiltonian (1)]. For such x another 
mechanism of the dip formation in the damping comes into effect. The 
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interaction constant 7k' — 7k'+k in Eq. (14) vanishes in the so-called "hot 
spots" - the crossing points of the Fermi surface with the boundary of the 
magnetic Brillouin zone. For small frequencies and k = Q the nearest 
vicinity of these points contributes to the damping. Consequently, the 
damping has a dip at Q. With the inclusion of the hole transfer to more 
distant coordination shells the interaction constant is changed, however 
the conclusion about the dip at the antiferromagnetic momentum in the 
damping remains unchanged. 

The dip disappears when the hole damping rj exceeds the frequency 
co > 0. Since the main contribution to the damping is made by states 
with energies — uj < ev — fi < and < £k'+k — /•* < w, this means that 
the above consideration is valid when there exist well defined quasipar- 
ticle excitations near the Fermi surface. 

Results of our calculations in a 120x120 lattice and experimental 
data are shown in Fig. 5. In agreement with experiment [20, 22] in 



0.08-, 




0.3 0.4 0.5 0.6 0.7 0.8 0.0S 0.10 0.15 

k(r.lu.) x 



Figure 5. (a) The dynamic structure factor vs. wave vector along the edge of the 
Brillouin zone for x = 0.015 (solid line), 0.06 (dashed line), and 0.12 (short dash). 
For all curves T = 58 K and lo — 2 meV. Inset: the structure factor measured 
[20, 21] in La 2 - a; Sr I Cu04 for x = 0.03, T = 2 K, lo = (filled squares and solid 
line), £ = 0.06, T = 25 K, oj = 2 meV (triangles and dashed line), and x = 0.12, 
T — 31 K, lo — 2 meV (open squares and short dash). Fitting curves in the inset are 
from Refs. [20, 21]. (b) The incommensurability parameter S vs. x for T — 58 K and 
lo — 2 meV. Inset: experimental data [20] for La2- a; Sr I Cu04. The wave vector and 
the incommensurability parameter are given in the reciprocal lattice units 2-k /a. 

our calculations peaks in the dynamic structure factor S(kuj) = [1 — 
exp(— uj /T)}~ l x" (koj) along the edge of the Brillouin zone [the direction 
(it, 7r) — (0, it)] are more intensive than those along the diagonal [the di- 
rection (it, it) — (0,0)]. As seen from Fig. 5, the calculations reproduce 
correctly the main known features of the incommensurate response ex- 
cept that the experimental 5 is somewhat larger than the calculated one. 
In experiment with increasing T the incommensurability decreases and 
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finally disappears. The calculations reproduce this peculiarity also. In 
the calculations the concentration dependence of 5 is mainly determined 
by the spin gap parameter A in uj^ in Eq. (9). For x < 0.12 it grows 
with x and saturates at larger concentrations [7] . 

With increasing uj the dip in the spin excitation damping is shallowed 
and finally disappears. Besides, on approaching loq the denominator 
in Eq. (10) will favor the appearance of the commensurate peak at Q. 
Thus, the low-frequency incommensurate maxima of x" converge to the 
commensurate peak at u « ujq. The dispersion of the maxima in x" 
above cjq is determined by the denominator in Eq. (10) and is close 
to that shown in Fig. 1. Consequently, the dispersion of the maxima 
in x" resembles two parabolas converging near the point (Q,cjq). The 
upper parabola with branches pointed up reflects the dispersion of spin 
excitations, while the lower parabola with branches pointed down stems 
from the momentum dependence of the spin excitation damping. Such 
kind of the dispersion is indeed observed in cuprates [3, 23]. 

Generally the incommensurate magnetic response is not accompanied 
by an inhomogeneity of the carrier density. In works based on the stripe 
mechanism [23] the magnetic incommensurability is connected with the 
appearance of a charge density wave. Notice however, that the magnetic 
incommensurability is observed in all lanthanum cuprates, while the 
charge density wave is detected in neutron scattering in those cuprates 
which are in the low-temperature tetragonal or less-orthorhombic phases 
(La2- a; Ba a ;Cu04, L^-y-rrNdySr^CuO^) [24]. In these phases the charge 
density wave is stabilized by the corrugated pattern of the in-plane lattice 
potential. It can be supposed that the magnetic incommensurability may 
be a precursor rather than a consequence of the charge density wave. 

6. Concluding remarks 

In this paper the projection operator technique was used for inves- 
tigating the magnetic properties of the t-J model of cuprate high-T c 
superconductors. It was demonstrated that the calculations reproduce 
correctly the frequency and momentum dependencies of the experimen- 
tal magnetic susceptibility and its variation with doping and temper- 
ature in the normal and superconducting states in YBa2Cu3C , 7_j / and 
lanthanum cuprates. This comparison with experiment allowed us to 
associate the resonance peak in YBa2Cu3C , 7_j / with the magnon branch 
modified in the short-range antiferromagnetic order. The lack of the res- 
onance peak in L^-ajSr^CuO^ was connected with an increased damping 
of spin excitations at the antiferromagnetic wave vector Q = (ir, tt). One 
of the possible reasons for the overdamped excitations is an increased 
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hole damping in this crystal. It was shown that for low frequencies the 
susceptibility is peaked at incommensurate momenta (tt ± 2tt5, tt) and 
(tt, tt±2tt5) . The incommensurability is the consequence of the dip in the 
momentum dependence of the spin excitation damping at Q. The dip ap- 
pears due to the Fermi surface nesting for low hole concentrations x and 
due to small hole-magnon interaction constants for moderate concentra- 
tions. In agreement with experiment for x < 0.12 the incommensurabil- 
ity grows nearly proportional to x and tends to saturation for x > 0.12. 
This is connected with the concentration dependence of the spin exci- 
tation frequency at Q. Also in agreement with experiment the incom- 
mensurability decreases with increasing temperature. The dispersion of 
the maxima in the susceptibility resembles two converging parabolas. 
The upper parabola with branches pointed up reflects the dispersion of 
spin excitations, while the lower parabola with branches pointed down 
stems from the momentum dependence of the spin excitation damping. 
In the considered mechanism the magnetic incommensurability is not 
accompanied by the inhomogeneity of the carrier density. 

Acknowledgments 

This work was supported by the ESF grant No. 5548. 

References 

[1] P. Bourges, in The Gap Symmetry and Fluctuations in High Temperature Su- 
perconductors, edited by J. Bok, G. Deutscher, D. Pavuna, and S. A. Wolf 
(Plenum Press, 1998), p. 349; H. He, Y. Sidis, P. Bourges, G. D. Gu, A. Ivanov, 
N. Koshizuka, B. Liang, C. T. Lin, L. P. Regnault, E. Schoenherr, and B. Keimer, 
Phys. Rev. Lett. 86, 1610 (2001). 

[2] G. Aeppli, T. E. Mason, S. M. Hayden, H. A. Mook, and J. Kulda, Science 279, 
1432 (1997). 

[3] M. Arai, T. Nishijima, Y. Endoh, T. Egami, S. Tajima, K. Tamimoto, Y. Shio- 
hara, M. Takahashi, A. Garrett, and S. M. Bennington, Phys. Rev. Lett. 83, 608 
(1999); D. Reznik, P. Bourges, L. Pintschovius, Y. Endoh, Y. Sidis, T. Matsui, 
and S. Tajima, cond-mat/0307591. 

[4] F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 (1988). 

[5] H. Mori, Progr. Theor. Phys. 34, 399 (1965); A. V. Sherman, J. Phys. A 20, 
569 (1987). 

[6] J. H. Jefferson, H. Eskes, L. F. Feiner, Phys. Rev. B 45, 7959 (1992); A. V. Sher- 
man, Phys. Rev. B 47, 11521 (1993). 

[7] A. Sherman and M. Schreiber, Phys. Rev. B 65, 134520 (2002); 68, 094519 
(2003); Eur. Phys. J. B 32, 203 (2003). 

[8] J. Kondo and K. Yamaji, Progr. Theor. Phys. 47, 807 (1972); H. Shimahara 
and S. Takada, J. Phys. Soc. Jpn. 61, 989 (1992); S. Winterfeldt and D. Ihle, 
Phys. Rev. B 58, 9402 (1998). 



14 

[9] A. K. McMahan, J. F. Annett, and R. M. Martin, Phys. Rev. B 42, 6268 (1990); 
V. A. Gavrichkov, S. G. Ovchinnikov, A. A. Borisov, and E. G. Goryachev, Zh. 
Eksp. Teor. Fiz. 118, 422 (2000) [JETP (Russia) 91, 369 (2000)]. 

[10] D. Z. Liu, Y. Zha, and K. Levin, Phys. Rev. Lett. 75, 4130 (1995); N. Bulut 
and D. J. Scalapino, Phys. Rev. B 53, 5149 (1996). 

[11] J. Bonca, P. Prelovsek, and I. Sega, Europhys. Lett. 10, 87 (1989). 

[12] D. Forster, Hydrodynamic Fluctuations, Broken Symmetry, and Correlation 
Functions (W. A. Benjamin, Inc., London, 1975). 

[13] J. L. Tallon, C. Bernhard, H. Shaked, R. L. Hitterman, and J. D. Jorgensen, 
Phys. Rev. B 51, 12911 (1995). 

[14] B. Keimer, N. Belk, R. G. Birgeneau, A. Cassanho, C. Y. Chen, M. Greven, 
M. A. Kastner, A. Aharony, Y. Endoh, R. W. Erwin and G. Shirane, Phys. Rev. 
B 46, 14034 (1992). 

[15] N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 (1966); P. C. Hohen- 
berg, Phys. Rev. 158, 383 (1967). 

[16] A. Sherman, cond-mat/0409379. 

[17] D. K. Morr and D. Pines, Phys. Rev. Lett. 81, 1086 (1998). 

[18] A. Sherman and M. Schreiber, Phys. Rev. B 69, 100505(R) (2004); A. Sherman, 
phys. status solidi (b) 241, 2097 (2004). 

[19] X. J. Zhou, T. Yoshida, D.-H. Lee, W. L. Yang, V. Brouet, F. Zhou, W. X. Ti, 
J. W. Xiong, Z. X. Zhao, T. Sasagawa, T. Kakeshita, H. Eisaki, S. Uchida, 
A. Fujimori, Z. Hussain, and Z.-X. Shen, cond-mat/0403181. 

[20] K. Yamada, C. H. Lee, K. Kurahashi, J. Wada, S. Wakimoto, S. Ueki, H. Kimura, 
Y. Endoh, S. Hosoya, G. Shirane, R. J. Birgeneau, M. Greven, M. A. Kastner, 
and Y. J. Kim, Phys. Rev. B 57, 6165 (1998). 

[21] S. Wakimoto, G. Shirane, Y. Endoh, K. Hirota, S. Ueki, K. Yamada, R. J. Bir- 
geneau, M. A. Kastner, Y. S. Lee, P. M. Gehring, S. H. Lee, Phys. Rev. B 60, 
R769 (1999). 

[22] H. Yoshizawa, S. Mitsuda, H. Kitazawa, and K. Katsumata, J. Phys. Soc. Jpn. 
57, 3686 (1988); R. J. Birgeneau, Y. Endoh, Y. Hidaka, K. Kakurai, M. A. Kast- 
ner, T. Murakami, G. Shirane, T. R. Thurston, and K. Yamada, Phys. Rev. B 
39, 2868 (1989). 

[23] J. M. Tranquada, H. Woo, T. G. Perring, H. Goka, G. D. Gu, G. Xu, M. Fujita, 
and K. Yamada, Nature 429, 534 (2004). 

[24] M. Fujita, H. Goka, K. Yamada, and M. Matsuda, Phys. Rev. Lett. 88, 167008 
(2002). 



